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The effect of a magnetic field on the dipole echo amplitude in glasses at temperatures of about 
10 mK caused by nonspherical nuclei with electric quadrupole moments has been studied theoreti- 
cally. It has been shown that in this case, the two-level systems (TLS's) that determine the glass 
properties at low temperatures are transformed into more complicated multilevel systems. These 
systems have new properties as compared to usual TLS's and, in particular, exhibit oscillations of 
electric dipole echo amplitude in magnetic field. A general formula that describes the echo amplitude 
in an arbitrary split TLS has been derived with perturbation theory. Detailed analytic and numer- 
ical analysis of the formula has been performed. The theory agrees qualitatively and quantitatively 
well with experimental data. 

PACS numbers: 61.43.-j, 76.60.Gv, 76.60.Lz, 81.05.Kf 



I. INTRODUCTION 

Glasses at temperature below 1 K display a number 
of universal properties which are fundamentally different 
from the properties of similar crystals. These properties 
are almost independent of the chemical composition and 
are therefore due to the structure of amorphous solid; 
more exactly due to absence of long-range order in the 
glass [l| . One of such universal phenomena is a two-pulse 
dipole echo. 

The essence of the effect is the following. When a glass 
is subjected to two short electromagnetic pulses at a fre- 
quency 1 GHz separated by a time interval r, a response 
at the same frequency may be detected at time r after 
the second pulse. 

According to the standard theory by Anderson, 
Halperin, Varma, and Phillips au these properties 
are associated with the existence of so-called two-level 
systems (TLS's) in the glass, which are the groups of 
atoms tuneling between two minima of a double-well po- 
tential. However, the microscopic nature of the TLS's in 
the majority of glasses remains unknown. 

An interesting phenomenon was discovered in 2002 Q . 
The amplitude of the dipole echo in some nonmagnetic 
glasses exhibited a strong nonmonotonic (oscillating) de- 
pendence on the magnetic field even in weak fields (about 
10 mT). Wurger, Fleischmann, and Enss Q suggested 
that the nonmonotonic magnetic field dependence of the 
dipole echo amplitude is caused by the presence of nuclei 
with the electric quadrupole moment (or, equivalcntly, 
the nuclear spin J > 1) in the glass. 

This suggestion qualitatively agrees with the experi- 
mental data, in particular, with the results of the recent 
measurements of the dipole echo in glycerol (C3H8O3) 
In these experiments, it was shown that the replace- 
ment of hydrogen, which has a nuclear spin J = 1/2 and, 
consequently, a zero quadrupole moment, by deuterium 



(J = 1 and a nonzero quadrupole moment) results in 
a more than an order of magnitude enhancement of the 
magnetic-field dependence of the echo amplitude. 

The magnetic-field dependence of the echo amplitude 
was theoretically studied by Wurger, Fleischmann, and 
Enss @, Q and by one of us However, a detailed 
comparison of the theoretical results with experiments 
was not carried out in those works. This work is aimed 
at the detailed comparison of the theoretical results 6] 
with the available experimental data. We perform a the- 
oretical analysis of the derived general formula in vari- 
ous limiting cases and show that the theory agrees with 
the experimental data reasonably well, both qualitatively 
and, in some cases, quantitatively. 



II. INTERACTION OF TWO-LEVEL SYSTEMS 
WITH NUCLEAR QUADRUPOLE MOMENTS 



Let us briefly recall how the presence of atoms with nu- 
clear quadrupole moments affects the properties of two- 
level systems and, in particular, the amplitude of the 
dipole echo in glasses @. Let one of the atoms that 
are displaced under the tunneling of the TLS have a nu- 
clear quadrupole moment. The energy of its quadrupole 
interaction with the microscopic electric field and Zee- 
man interaction with the external magnetic field leads to 
TLS level splitting, so that the TLS becomes a system 
with two identical level sets each containing 2 J+ 1 levels. 
These sets are separated by the energy E of the original 
TLS, which is much greater than the energy splitting 
within the set. 

The Hamiltonian of such a system may be written in 
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the form 
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Here, 1 CT and 1q are the 2 x 2 and TV x N (N = 2 J + 1) 
unity matrices in the spaces of the TLS and the pro- 
jections of the nuclear spin, respectively. A is the dif- 
ference between the double-well potential minima. Ao 
is the tunneling amplitude of the initial (unsplit) TLS, 
E = ^/A 2 + A 2 , is the total energy of the TLS disre- 
garding the fine structure effects. F is the ac electric 
field, and m is the electric dipole moment of the TLS. 

Wq = M • H + Q a pip a p(0), where M is the nuclear mag- 
netic moment, H is the external magnetic field, Q a /3 is 
the operator of the nuclear quadrupole moment [7[ (ex- 
pressed in terms of the nuclear spin operator), and <p a p(Q) 
is the second-derivative tensor of the electric potential 
taken at the double-well potential maximum. 

The second term in Eq. (JTJ), which contains Wq, causes 
the fine splitting of the TLS. We assume that the ba- 
sis nuclear spin wavefunctions are eigenfunctions of the 

operator Wq and its matrix is diagonal in this basis. 
The last term in the Hamiltonian ([T]) includes the oper- 
ator Vq = Q a 0<p a p(O)\x min \/x o , where (f' a/3 (0) — is the 
derivative of the tensor <p a p(x) with respect to the gener- 



alized coordinate x of the TLS at 



and 



IXq is 



the ratio of the atom displacement under the tunneling 
of the TLS to the characteristic interatomic distance xo. 

The first term in Eq. fl} determines the two-level sys- 
tem, the second one is responsible for the fine splitting of 
the TLS due to the magnetic and quadrupole moments 
of the nucleus, and the third term determines the matrix 
elements of the transitions between various levels in the 
microwave field. Finally, the fourth term couples the split 
levels of the system, which leads to allowed transitions 
between different fine structure levels and, ultimately, to 
the oscillations of the dipole echo in the magnetic filed. 
Since \x m i n \/xo <C 1, the last term is small and may be 
taken into account perturbatively. 

Performing the calculations in the lowest (second) or- 
der of perturbation theory, we come (similar to Q) to 
the following expression for the echo amplitude: 
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Here, r is the time interval between two excitation pulses, 
£nm = (W Q ) nn - (Wq) toto is the energy differences be- 
tween the fine structure levels of the TLS. Equation ([2|) 
differs from similar Eq. (44) in [6] basically due to of 

2 



the absolute-value square 



(Vq), 



which implies that 



our consideration is not limited by only the real-valued 
matrix elements (VQ) nm . 

The magnetic field enters Eq. ([2|) first via the energy 
difference e nm and, second, implicitly via the matrix el- 
ements iyo) n mi owing to the fact that they were calcu- 
lated in the basis in which Wq is diagonal. 

According to our numerical calculations, the latter de- 
pendence is usually insignificant for finding the echo am- 
plitude for the nuclei with integer spin J. However, when 
the spin is a half-integer and the levels n and m are de- 
generate according to the Kramers theorem, the matrix 
element (Vg)„ m cannot remove degeneracy and, conse- 
quently, must vanish in a zero magnetic field. 

Formula $2$ allows a numerical calculation of the echo 
amplitude for given parameters of the system. However, 
a direct application of this formula is complicated by the 
necessity of solving the algebraic equation of degree 2 J+l 
with subsequent averaging of the results (which cannot 
be done in the general form even for J = 1). Thus, 
Eq. ([2|) must be investigated in various limiting cases 
and analyzed numerically. 



III. ANALYSIS OF LIMITING CASES 

We begin the analysis of Eq. @ with three indepen- 
dent energy scales appearing in it. The first two are the 
Zeeman energy Eh and the energy Eq of the quadrupole 
interaction of the nucleus with the microscopic electric 
field. The third scale E T = H/t is determined by the 
time interval r between the pulses. 

When the magnetic field is so high that Eh is greater 
than the other two scales, we can replace sin 4 (e„ m T/2fi.) 
by its average value and neglect the quadrupole interac- 
tion energy in the calculation of Wq (and consequently 
Enm)- In this case, the echo amplitude turns out to be 
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(3) 



where the coefficient C is independent of the magnetic 
field. 

Thus, the echo amplitude asymptotically (as 1/H 2 ) 
approaches a constant value with an increase in the mag- 
netic field, which agrees with the experimental data for 
high fields [3, H| . A less trivial magnetic-field dependence 
of the echo amplitude takes place when the Zeeman en- 
ergy Eh is comparable with (or less than) at least one of 
the other energy scales Eq or E T . 
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A. Small quadrupole splitting 

Consider now the case where the quadrupole interac- 
tion energy is Eq <C Eh,E t . In this case, Eh and 
E T may be comparable. Then, we can neglect the 

quadrupole energy in the calculation of Wq (but we 

still have to retain it in Vq). Correspondingly, the 
set of the eigenfunctions of the projections of spin J 
on the direction of the magnetic field can be used as 
the basis. In this case, the fine structure levels are 
E n = /ii?(J — n — 1)/ J, where /i is the nuclear magnetic 
moment and n — 1,2, 2J + 1. 

In this case, since glass is isotropic and the trace of the 
tensor <p' a p(0) is zero, the average squares of the absolute 



values of all matrix elements 



(V Q )r 



can be expressed 



in terms of a single constant. Then, the dipole echo am- 
plitude (as a function of the magnetic field) is determined 
by the formula 
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Here, C is the coefficient independent of the magnetic 
field (and different from C in Eq. ©). 

The magnetic field enters into Eq. (U|) as the product 
fiHr. This implies that the variation of the dipole echo 
amplitude with H should scale as 1 /r along the horizon- 
tal axis. 




FIG. 1: Magnetic-field dependence of the dipole echo ampli- 
tude in the BK7 glass (in arbitrary units) for r = (from top 
to bottom) 0.75, 1.5, 2, 3, 6, and 12 /is. The points are the 
experimental data from || and solid lines are the theoreti- 
cal curves for the nuclear magnetic moment of 10 B calculated 
with the use of Eq. @ . 



Let us compare Eq. ([4} with the measurement of the 
echo amplitude in the BK7 glass as a function of the 
magnetic field H (Fig. [T]). Obviously, the experimental 
dependence of the echo amplitude on H has two compo- 
nents, one of which keeps its scale with a change in r 
and the other one changes the horizontal scale as 1/t. 
The second component is well described by Eq. (U) with 
/i equal to the nuclear magnetic moment of the 10 B atom 
(which is present in the BK7 glass and has a nuclear spin 
J = 3). We may suggest that the nonscalable compo- 
nent of the dependence is determined by the contribu- 
tion of other atoms (e.g., Na or 11 B) with the nuclear 
quadrupole moment, for which the condition E T 3> Eq 
is not fulfilled. 



B. Magnetic field dependence of echo amplitude in 
weak fields 



Consider now another limiting case where Eq ~ E T 
but the Zeeman energy is low, Eh *C Eq,E t . This will 
give us an idea of the behavior of an echo amplitude in 
small fields. 

At H — 0, the energy differences e nm between the fine 
structure levels of the TLS are fully determined by the 
energy of the quadrupole interaction of the nucleus with 
the internal field and enter the final formula for the echo 
amplitude through the function 
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The plot of this function (see Fig. [21 left panes) is a set 
of peaks with the l/y 2 envelope. The character of the 
magnetic-field dependence of the dipole echo amplitude 
in weak fields is determined by the fact that the values 
of yo corresponding to the zero-field splitting may fall to 
a minimum or maximum of dependence ([5]). 

Consider for example the simplest case of the nucleus 
with J = 3/2. This nucleus has four fine structure levels 
and, respectively, three independent energy differences 
E nm . However, two of them (e\2 and £34) vanish in zero 
magnetic field, according to the Kramers theorem. The 
corresponding matrix elements (Vo) nm also vanish. As 
a result, the contribution of these two pairs of levels ap- 
pears to be small compared to the contribution of the 
remaining e 2 3 pair. 

Thus, in this case, the behavior of an echo in weak 
fields is determined by the single energy difference £23. 
Suppose that the time interval r between the pulses is 
such that yo = £23 (H = 0)r/2h — irk with integer fc, 
which corresponds to a minimum of function ((5|). Then, 
in the presence of the magnetic field, y must shift from 
the minimum position and, according to Eq. ([5]), the echo 
amplitude must decrease, irrespective of the properties of 
the system. Therefore, the magnetic-field dependence of 
the dipole echo should have a maximum at H = 0. 

Figures [^b and show the numerical results for the 
echo amplitude calculated with the use of Eq. ^ for the 



4 



a) 0.6 




FIG. 2: Positions of yo and the respective results of the sim- 
ulation of the dipole echo amplitude for J = 3/2. 



case of J = 3/2 and yo = ^ and 2w, respectively (contri- 
butions of all levels were taken into account). Clearly, in 
agreement with the prediction, the echo amplitude has a 
maximum at zero magnetic field. 

On the contrary, if r is such that yo corresponds to a 
maximum of function 1(5]). then the echo amplitude must 
increase in the presence of the magnetic field and the 
magnetic-field dependence of an echo should have a min- 
imum at H = 0. The corresponding numerical results 
are shown in Figs. [2k and[2fc. 

For an arbitrary spin, this analysis is possible if the 
contribution of one energy splitting e nm prevails over 
the contributions of the other pairs of levels. This oc- 
curs, e.g., if the other s nm values are large and can be 
disregarded in Eq. ^ due to the term l/e^ m . 

Thus, it turns out that a zero-field maximum of the 
echo amplitude should change to a minimum and vice 
versa with a change in the time interval r between the 
pulses. 
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FIG. 3: Experimental magnetic-field dependence of the dipole 
echo amplitude in glycerol 0, Qi| with r = (a) 3.5, (b) 7, and 
(c) 14 /is and (d) the respective zero-field positions y a , and 
Vc of yo- 



A similar behavior was experimentally observed in 
glycerol with hydrogen replaced by deuterium [To| . 
Deuterium has spin 1, which corresponds to three ze- 
rofield energy splittings £12 oc 2r), £23 oc 3i] and £13 oc 
3 + 77, where < r\ < 1 is the parameter of the micro- 
scopic field asymmetry. If r\ <C 1, the first level splitting 
is much smaller than the other two and, in agreement 
with the above analysis, may play the main role in the 
magnetic-field dependence of the echo amplitude. 

Figure 3 presents the experimental results [H, [l(| and 
the positions of yo for different values of r. They were 
calculated with the fitted value en/h « 110 kHz (which 
is comparable with an experimental value of 150 kHz [4J]). 
Clearly, the echo amplitude as a function of the magnetic 
field has a minimum at H = when yo appears near the 
maximum of function ([5]) and vice versa; i.e., the echo 
amplitude has a maximum at zero field when yo appears 
near the minimum of the function |(5J). 

Summarizing, the elaborated theory of the dipole echo 
in the magnetic field with the inclusion of the nu- 
clear quadrupole moments provides the explanation of 
all characteristic features observed in experiments and, 
in some cases, yield even a quantitative agreement be- 
tween the experiments and analytic calculations. How- 
ever, the comparison of the magnetic-field dependence 
of the dipole echo amplitude in an arbitrary glass with 
the theory is complicated by superimposed contribu- 
tions from the atoms with different nuclear magnetic and 
quadrupole moments. 
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